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SECTION  I 


INTRODUCTION 


Suppression  of  sympathetic  detonation  between  stored  munitions  has  been  a  very  important  issue  in 
the  1990  s  for  all  branches  of  the  Department  of  Defense. 

During  the  Air  Force's  Insensitive  High  Explosives  (IHE)  development  program,  a  series  of  live  typical 
munition  storage  items  in  a  standard  storage  configuration  was  tested  using  AFX- 1 100,  TNTAvax/AI 
(66/16/18  percent  by  weight).  The  storage  configuration  is  shown  in  Figure  1.  For  symmetry  and  worst 
case  conditions,  the  donor  was  placed  in  the  bottom  middle  position.  It  was  found  that  the  left  and  right 
bottom  and  top  center  items  did  not  detonate  when  exposed  to  the  detonation  of  the  donor.  It  was  also 
observed  that  the  left  and  right  diagonal  items  consistently  detonated.  Since  the  items  did  not  detonate  in 
a  side  by  side  test  at  the  same  diagonal  distance,  it  was  hypothesized  that  the  confinement  of  the  donor 
item  due  to  the  top  center,  bottom  left  and  right  items  was  causing  an  enhancement  of  the  munition  case 
velocity  up  to  the  critical  initiation  pressure  of  AFX-!  100  for  items  located  in  the  diagonal  position. 


Figure  I.  Typical  Munition  Storage  Configuration 

A  second  series  of  tests  was  conducted  to  verify  the  hypothesis.  The  tests  were  designed  to  alleviate 
some  of  the  confinement  of  the  donor  by  elevating  the  top  row  of  munition  items.  The  minimum 
separation  distance  for  the  munition  case  was  13  mm  horizontally  and  vertically.  The  diagonal  distance 
was  133  mm.  As  the  top  row  of  items  were  moved  up,  which  reduced  the  confinement,  the  diagonal  item 
no  longer  sustained  a  detonation.  It  appeared  that  the  hypothesis  was  uue.  But  upon  further  investigation 
using  hydrocodes  (Reference  1)  it  was  determined  that  the  casewall  of  the  donor  item,  with  the  top  row  of 
items  elevated,  produced  a  higher  velocity  during  the  detonation  than  did  the  standard  confinement  case. 
Figure  2  is  a  hull  calculation  1(X)  ps  into  the  detonation  of  the  donor  with  the  items  in  the  normal  stacking 
configuration.  Figures  3  and  4  are  time  history  data  of  donor  casewall  velocity  and  acceptor  pressure  for 
this  calculation. 
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Figure  2.  Hull  Calculation  No.  1 


Figure  3.  Hull  Calculation  No.  I  Showing  Casewall  Velocity  at  Impact  of  the  Donor  Casewall 

Notice  that  the  flyer  casewall  velocity  for  the  suindord  configuration  case  was  1.55  mm/us  and  the 
pressure  induced  inside  the  acceptor  explosive  was  55  kbars.  By  way  of  comparison,  the  critical  initiation 
pressure  for  AFX-1 100  as  measured  by  the  modified  Expanded  Large-Scale  Gap  Test  (ELSGT)  is  between 
53  and  56  kbars.  The  ELSGT  pulse  duration  is  very  similar  to  that  calculated  for  the  diagonal  item  shown 
in  Figure  2.  Thus  because  the  calculation  predicts  that  the  diagonal  item  is  at  the  initiation  threshold  for 
AFX-1 100,  we  should  expect  the  acceptor  to  detonate. 

2 


Impact 


The  next  series  of  calculations  were  performed  at  a  non-detonating  height  for  the  diagonal  item  of  76 
mm  as  measured  vertically  on  the  outside  of  the  items.  Notice  in  Figure  S  that  at  100  (is  the  flat  plate 
generated  from  the  donor  casewall  appears  to  have  thinned  more  than  in  the  previous  test  (see  Figure  2). 
Thinning  of  the  casewall  is  directly  related  to  the  amount  of  expansion  the  munition  case  is  allowed  to 
undergo.  As  a  general  rule  it  is  assumed  that  the  munition  case  will  expand  up  to  two  times  the  initial 
radius  before  it  breaks  up.  In  Figure  6  the  pressure  induced  inside  the  acceptor  item  is  calculated  to  be 
44,6  khars,  approximately  10  kbars  below  critical  initiation  pressure.  The  velocity  of  the  thinned  casewall 
at  impi.  '  an  the  diagonal  acceptor  item  is  1.62  mm^ts  as  shown  in  Figure  7.  The  calculation  predicts  no 
initiation  in  this  instance  and  is  consistent  with  the  experimental  observation. 


Figure  5.  Hull  Calculation  No.  2 
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Figure  6.  Hull  Calculation  No.2  Showing  Pressure  Pulse  Inside  of  Acceptor  Item 
Due  to  Donor  Impact 

As  the  top  row  of  bombs  is  raised,  the  donor  casewall  expands  further,  hence  thins  more,  prior  to 
impact  with  the  center  item  of  the  top  row.  Although  confinement  docs  not  result  in  higher  case  velocity, 
it  does  reduce  case  thinning  given  the  case  is  not  free  to  expand  radially. 


Figure  7.  Hull  Calculation  No.  2  Showing  Casewall  Velocity  at  Impact  of  the  Donor 
Casewall 

A  calculation  was  ai  c  performed  to  determine  if  the  complex  geomeuy  depicted  in  Figure  2  could  be 
modeled  by  impacting  a  single  cylinder  with  a  fiat  plate  as  shown  in  Figure  8.  The  pressure  ob.scrvcd  at 
the  first  unmixed  cell  was  55.7  kbar;  recall  that  the  pressure  predicted  for  condition  in  Figure  2  was  55 
kbar. 
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Figure  8.  Hull  Calculation  of  Flyer  Plate  With  Pressure  Pulse  Signal  Induced  in  the 
Acceptor 

The  next  calculation  was  performed  to  see  if  the  detonation  producLs  contributed  to  the  ove  ^il  energy 
of  the  Hyer  plate.  We  have  hypothesized  that  the  munition  test  as  shown  in  Figure  1  is  a  long  impulse 
event.  However,  from  Figure  4,  very  little  area  exists  under  the  initial  pressure  pulse.  This  implies  that 
the  pressure  duration  is  coniiollcd  by  the  thickness  of  the  impacting  cascw  alls  with  little  contribution  from 
the  detonation  products.  Based  on  the  calculation  of  the  detonating  donor  item,  at  impact,  the  gases  have 
expanded  into  a  volume  VA'o  of  between  2  and  3.  A  complete  history  of  the  expansion  isentrope  of  the 
donor  item  is  shown  in  Figure  lU.  The  pressure  associated  with  this  expansion  is  between  2  and  5  kbars. 
To  verify  that  the  pressure  behind  the  flyer  plate  was  not  contributing  to  the  input  pressure  as  seen  by  the 
acceptor  explosive,  a  calculation  shown  in  Figure  9  was  performed  with  10  kbars  of  pressure  behind  the 
flyer  plate.  All  the  other  conditions  were  kept  the  same.  The  pressure  at  the  first  unmixed  cell  inside  the 
acceptor  case  was  58  kbars.  It  is  apparent  that  the  contribution  of  the  product  gases  lo  the  input  pressure  is 
minimal.  Based  upon  all  of  the  information,  it  is  believed  the  diagonal  acceptor  is  being  initiated  by  a 
shock  to  detonation  transition  that  is  directly  related  to  the  velocity,  thickness,  and  shape  of  the  flyer  plate. 
From  this  information  the  Explosively  Driven  Flyer  Plate  Test  was  developed  to  reproduce  the  boundary 
conditions  as  observed  with  experiment  and  hydrocode  calculations. 


X  (cm)  TVne  (us) 


Figure  9.  Hull  Calculation  of  the  Detonation  Product  Gas  in  Conjunction  With  the  Flyer  Plate 
and  the  Pressure  Pulse  Calculation  for  the  Inside  of  the  Acceptor  Bomb 
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SECTION  U 


SIMPLIFIED  MODEL  DEVELOPMENT 

The  development  of  a  simplified  model  to  predict  the  results  of  the  flyer  plate  test  was  considered 
essential  for  numerous  reasons.  Principally,  if  the  event  was  well  understood,  a  model  could  be 
developed  and  also,  with  a  simplified  model,  the  necessity  of  understanding  both  the  equations  of  state 
and  numerical  methods  employed  in  a  large,  rnultimaterial.  multidimensional  finite  difference  code,  such 
as  the  OTI*HULL  code,  could  be  avoided.  This  is  particularly  true  of  resolution  requirements  that  in  a 
simplified  model  can  be  avoided  altogether.  For  the  development  and  checkout  of  the  model,  the 
OTI*HULL  code  was  used. 

Tiie  model  is  broken  into  three  pans:  (1)  the  cylinder  expansion,  (2)  impact,  casewall  pressure 
divergence,  and  (3)  the  acceptor  explosive  evaluation.  Each  of  these  parts  will  be  addressed  separately. 

l.Cyi.INDER  EXPANSION 

The  model  currently  allows  the  user  three  separate  means  of  determining  the  donor  casewall  impact 
velocity.  The  first  method  is  a  straightforward  integration  of  the  expansion  isentrope.  The  second 
method  attempts  to  better  match  cylinder  expansion  data  through  the  use  of  impulse  momentum.  The 
last  is  a  combination  of  the  first  two  methods.  The  Jones,  Willons,  and  Lee  (JWL)  descr^on  of  the 
cqiansion  isentrope  is 

=  +  +  (1) 
where  a,b,c,rj,  and  r2  are  constants  which  allow  fora  non-constant  y  and  V  which  is  the  ratio  of  VA'o 

The  integration  of  the  expansion  isenuope  can  be  found  in  numerous  references  but  for  this 
application  the  method  given  by  Miller  (Reference  2}  in  Equation  1  is  most  suitable.  To  avoid  repetition, 
only  a  summary  will  be  given  here.  The  energy  as  a  result  of  Pdv  work  is  written  as: 


E  V 

J  =  J  pdV  (2) 

£«  1 

or 


E  =  E(\/  =  1)  -  +  ~  V-) .  (3) 

n  ri  w 

From  the  Gurney  relation, 

£^=:E„(y)(l  +  0.5p/y)  (4) 

v  =  yl2Eg). 
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Then  setting  E(V=1)  equal  to  the  heat  of  detonation,  the  energy  given  by  Equation  3  can  then  be 
substituted  into  Equation  4  for  Ejg  giving  the  velocity  as  a  function  pC  volume  expansion. 

A  plot  using  this  method  can  be  seen  in  Figure  1 1  with  a  hydrocodc  calculation  also  shown  for 
comparison.  Late  time  data  is  fit  very  well  but  early  time  response,  which  is  controlled  by  momentum 
Uansfer,  is  not. 

I  WCM  CO  CntNOtB 


To  overcome  this  difficulty,  a  model  based  upon  impulse  and  momentum  as  given  by  Equation  5  was 
developed. 

Mdv  (5) 


where  F  =  force 

dt  =  loading  duration 
M  =  mass 

dv  s  change  in  velocity. 

Substituting  the  pressure,  inside  casewall  area  and  mass  into  Equation  5  yields; 


V,  =  Vi - 


.  FA.-,dt 

^  ■ 

M 


(6) 


Aj.|  ®  Area  of  the  inside  ca.scwall 

Since  this  equation  by  itself  does  not  yield  a  unique  result,  all  quantities  are  sealed  to  the  1-inch 
copper  cylinder  test.  Tb  '  equations  arc  cast  in  plane  geometry  so  that  dt  is  taken  as  the  transit  time  of  the 
detonation  wave  over  1  cm  00  mm).  This  result  is  further  broken  down  into  n  divisions.  Further,  the 
initial  casewall  volume  (V/Vq  =  1)  which  implies  that  the  Chapman  Jouget  (CJ)  pressure  is  never 
attained.  To  avoid  this  difficulty,  the  CJ  volume  is  calculated  as  Equation  7; 
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(7) 


Vy  = 


T 

Y  +  1 


Vo. 


This  value  is  subtracted  from  the  initial  volume  ratio  (1)  so  that  the  volume  ratio  input  to  the  JWL 
equation  of  state  (EOS)  is  reduced  for  a  time  by  this  amount.  However,  once  the  reaction  zone  passes  a 
given  part  of  the  cylinder,  this  volume  reduction  no  longer  applies.  It  is  assumed  that  when  the  shock 
wave  reaches  the  cylinder  free  surface,  this  volume  reduction  is  removed.  To  calculate  the  shock 
velocity  in  the  casewall  due  to  the  impact  of  a  detonation  wave,  Equation  8  is  applied. 


u,~u  f  2 

y^/(y  +  1)J  yiy  +  i)  p/ 

where  u  i ,  P  i ,  D,  7  refer  to  the  CJ  condition. 

u  =  particle  velocity.  P  =  pressure,  D  =  detonation  velocity. 

The  intersection  of  Equation  8  with  the  casewall  Hugoniot  provides  the  interface  particle  velocity. 
Substitution  of  the  interface  particle  velocity  into  the  casewall  Hugoniot  provides  the  shock  velocity  in 
the  casewall. 

Knowing  the  time  step  (divided  by  n  divisions),  the  number  of  required  iterations  can  be  calculated 
as; 

Iter  =  (cwi/US(;)(n/dcji)  (9) 

where 

cwi  s  casewall  thickness 
USg  =  shock  velocity  in  casewall 
n  =  number  of  divisions 
deji  =  detonation  velocity  of  explosive. 


Since  the  equations  are  cast  in  plane  geometry,  only  radial  expansion  can  be  considered.  However,  it 
has  been  found  necessary  to  allow  the  volume  to  expand  axially  as  a  function  of  the  panicle  velocity. 
This  can  be  expressed  as: 

h  =  hi.i  +  Up*k  (10) 

The  particle  velocity  (Up)  can  be  calculated  from  the  momentum  equation  as: 

Up  =  P/PoUs  (11) 

where  is  the  ambient  casewall  density 

The  shock  velocity,  Us,  is  given  as  Equation  12. 
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The  variable  Y  can  be  calculaicd  from  (he  JWL  EOS  as: 


Y  =  v(rjae'^''^  +  r2bc*^'^+  c(w+l)vC'''+2) 

-  .  (13) 


To  match  the  cylinder  expansion  with  a  minimum  number  of  time  steps,  the  interface  pressure  may 
also  be  applied  for  a  short  duration  but  not  exceeding  the  number  of  iterations  given  by  Equation  10. 


With  these  equations,  the  casewall  expansion,  hence  volume  increase,  can  be  calculated  from 
Equation  6.  The  number  of  divisions  (n)  and  the  axial  expansion  constant  k  are  determined  from  1-inch 
copper  cylinder  test  results  (or  hydrocode  calculations).  Both  constants  are  varied  until  (he  calculated 
cylinder  test  matches  the  experimental  results.  Values  for  n  vary  between  10  and  30  and  is  essentially 
controlled  by  the  explosive  energy.  Generally,  the  more  energetic  the  explosive  the  more  divisions 
required.  The  axial  expansion  constant  is  usually  between  0.1  and  0.2.  Calculations  performed  at  other 
charge  to  mass  ratios  arc  scaled  to  the  1-inch  cylinder  results.  The  axial  expansion  constant  is  given  as: 


where 


k  =  ( 1  -  cm,  +  xp )( r^ )  ^14) 

cm^  =  ratio  of  charge  to  mass  of  a  1 -inch  copper  cylinder  of  the  same 

explosive  to  the  ratio  of  the  current  casewall  material  to  explosive 

xp  s  axial  expansion  constant  for  one  copper  cylinder 

tf  =  ratio  of  current  inside  radius  to  inside  radius  of  a  one  inch  copper 
cylinder. 


Note  (hat  this  term  reduces  to  zero  for  the  1-inch  test  with  no  axial  expansion  (xp).  Values  for  the 
axial  expansion  constant  arc  generally  in  the  range  of  0.1  to  0.2  for  the  1  -inch  test.  Results  for  this  model 
are  shown  in  Figures  12  and  13  with  the  first  being  a  copper  cylinder  test  and  the  second  a  steel  cylinder 
with  a  cross  section  corresponding  to  e  general  purpose  bomb. 


I  wch  cu  cniMOCR 


Figure  12,  1 -inch  Copper  Cylinder  Test  for  AFX- II 00 
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Figure  13.  Sled  Cylinder  Casewall  Velocity  Test  Using  APX-llOO 

Velocities  are  limited  to  the  Gurney  velocity.  At  a  volume  ratio  of  approximately  2,  the  impulse- 
momentum  calculation  can  be  terminated  and  the  remaining  expansion  calculated  ns  a  function  of  energy. 
This  method  is  preferred  since  the  energy  method  is  significantly  faster  and  produces  better  late  time 
results.  The  casewall  motion  as  a  function  of  displacement  is  now  determined.  The  casewall  thickness  as 
a  function  of  expanding  radius  can  be  calculated  through  the  continuity  equation  as; 

t  -t-4rF^^r77^  (15) 

CM 

where  r  =  the  current  radius 

fj  =  initial  inside  radius 
fg  «  initial  outside  radius. 

2.  IMPACT  PRESSURE  DIVERGENCE 

The  results  of  the  previous  section  indicate  the  impact  velocity  can  be  calculated  to  within  a  few 
percent.  From  the  impact  velocity,  the  initial  pressure  in  the  casewall  can  be  determined.  Without 
divergence,  that  pressure  will  be  uniform  throughout  the  material  thickness.  However,  through 
hydrocode  analysis,  a  significant  amount  of  pressure  divergence  can  be  observed  in  the  casewall.  Figure 
14  depicts  the  results  of  three  impact  scenarios  involving  a  purely  one-dimensional  flat  plate  impact,  a 
flat  plate  impacting  a  right  circular  cylinder,  and  a  curved  plate  impacting  a  right  circular  cylinder. 


Figure  14.  Three  Impact  Scenarios  Using  Hytlrocodcs 
II 


The  pressure  pulse  response  in  the  casewall  and  in  the  explosive  can  be  seen  in  Figure  IS.  Note  that 
the  pujse  duration  is  the  same  for  all  three  cases. 


o  t'Ttet  Of  fi>it  euw*ni»t  «  cascwmi 


Figure  15.  Pressure  Pulse  Response  in  the  Casewall  and  the  Explosive 

These  calculations  were  performed  using  the  lagrange  module  in  OTI*HULL  at  a  resolution  required 
to  predict  the  analytical  one-dimensional  flat  plate  pressure.  The  ratio  of  impacting  plate  thickness  to 
acceptor  plate  is  O.S.  The  impacting  velocity  is  l.S  mm/usec.  As  can  be  seen,  each  case  begins  at  the 
same  initial  pressure,  which  can  be  calculated  through  the  momentum  equation  and  material  Hugoniot. 
The  one-dimensional  momentum  equation  is  written  as: 

PapoUsUp  (16) 

where  po  is  the  initial  density 

Us  is  the  shock  velocity 
Up  is  the  panicle  velocity. 

For  most  solid  materials,  the  shock  v-  <iiy  can  be  related  to  the  particle  velocity  at  the  linearly  as: 

Us  =  C  +  SUp  (17) 

where  C  is  the  ambient  longitudinal  sound  speed 

S  is  the  Hugoniot  slope. 

For  like  materials,  the  interface  pressure  at  impact  can  be  expressed  as: 

P  =  (po/2)(Up/2)(C  +  S  Up/2).  (18) 

For  dissimilar  materials,  a  reflected  Hugoniot  Is  employed.  The  particle  velocity  at  the  interface  in 
the  receiving  material  can  be  expressed  as; 

Up=2U,-U.  (19) 

The  unsubscripted  variables  represent  the  interface  conditions.  Equations  16  and  19  can  then  be 
written  on  both  sides  of  the  interface  and  the  resulting  interface  particle  velocity  solved  for,  which  in  turn 
provides  the  Interface  pressure.  Figure  16  depicts  two  impact  scenarios  involving  different  plate 
thicknesses  and  impact  velocities. 
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Figure  16.  Lagrangian  Calculations  of  Flyer  Plate  Impacts 

Note  the  radius  of  curvature  of  the  resulting  pressure  wave  due  to  the  impact. 
begins  at  the  impact  point  and  terminates  at  the  end  of  the  impacting  plate.  The  magnitude  of  the  iniual 
shock  wave  pulse  can  be  reduced  in  two  ways;  either  through  wave  divergence  or  rarefactions  resulting 
from  free  surface  interaction.  The  problem  then  is  to  determine  the  minimurn  diameter  equivalent  plate 
that  reproduces  the  conditions  involved  in  the  original  impact  problem.  Work  involving  spherical 
fragments  (Reference  3)  indicated  that  the  equivalent  flat  plate  was  controlled  by  lateral  free 
relfef  that  occurred  when  the  plate  closing  rate  was  exceeded  by  the  velocity  of  rarefaction  waves,  whic 
were  set  equal  to  the  ambient  sound  speed  in  the  material.  Since  the  spheres  were 
dimensional  region  was  considered  to  exist  in  the  impacted  material  prior  to  relief  due  to  ^^facuo" 
waves.  For  a  sphere  impacting  a  plane  surface,  as  shown  in  Figure  17.  with  umc  referenced  to  the  cent 
of  the  sphere,  the  equivalent  plate  diameter  was  given  as: 


di  =  2y(l+Co2/uj2)J 


(20) 


X.ut 


Figure  17.  Sphere  Impacting  a  Plane  Surface 

The  geometry  of  the  event  is  shown  in  Figure  18.  Victor  (Reference  4)  extended  this  to  curved  plates 


with  the  equivalent  plate  given  as: 


where 


dj  =  2ra(v/Coc)(rd  +  L)/(ra  +  rj  +L) 
r_  is  the  radius  of  the  acceptor 


(21) 


Vf  is  the  velocity  at  donor  radius  r 
rjj  is  the  initial  donor  radius 

L  is  the  shortest  distance  between  donor  and  acceptor. 
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Work  performed  by  Orcen  (Reference  5)  suggested  that  lateral  rarefactions  enter  into  the  onc- 
dimensional  region  from  45  degree  angles.  Therefore  if  the  lateral  extent  was  known,  propensity  toward 
sympathetic  detonation  could  be  determined  by  examining  the  depth  into  the  explosive  the  one¬ 
dimensional  region  extended  and  comparing  that  value  to  the  explosive  critical  diameter.  For  depths  iess 
than  the  critical  diameter,  an  expanding  wave  approximation  of  the  form  (Equation  22)  could  be  used. 


where 


Peqv  =  Po 


1 

.  do  J 


xi  is  the  distance  into  the  acceptor 
do  is  the  equivalent  flat  plate  diameter. 


(22) 


Figure  18.  Geometry  Used  to  Describe  Equation  19 

As  mentioned  previously,  considerable  divergence  exists  within  the  cascwall  and  must  be  taken  into 
account.  In  addition,  rarefactions  also  enter  due  to  rear  surface  relief  and  may  conuol  the  plate  diameter. 
Hydrocodc  results  suggests  that  a  better  approximotion  to  the  equivalent  flat  plate  can  be  determined  by 
replacing  the  ambient  sound  speed  in  either  Equations  20  or  21  with  the  shock  velocity  appropriate  to  the 
impact  velocity  determined  through  the  momentum  equation.  This  result  is  then  compared  to  the  plate 
diameter  calculated  assuming  rear  surface  relief  is  the  controlling  mechanism.  An  initial  approximation 
of  the  available  time  is  made  by  determining  the  uansit  time  of  a  shock  wave  across  the  acceptor 
thickness.  This  time  is  compared  to  the  time  of  a  rearward  propagating  shock  wave  into  the  donor  plate. 
The  donor  plate  rear  surface  continues  at  the  initial  impact  velocity  until  encountering  the  rearward 
propagating  shock  wave. 
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This  determines  both  the  plate  thickness  and  pulse  duration  into  the  explosive.  This  time  can  be 
expressed  by; 

t  =  plate  thickness  /  ( impact  velocity  +  shock  velocity  )  (23) 

The  surface  then  moves 


rsm  =  t*v 


(24) 


where  v  =  impact  velocity,  and  rsm  =  rear  surface  movement. 

The  new  plate  thickness  is  then 

Tpltn  =  Tp  -  rsm  (25) 

where  the  old  plate  thickness  Tp  is  calculated  through  the  continuity  equation  as; 


tp-roi-  JTZ  -  (r»  -ri.) 


where  ro,  is  the  outside  radius  at  impact 

ro<i  is  the  initial  outside  radius 
ri,  is  the  initial  inside  radius . 


Initial  efforts  to  determine  the  extent  of  pressure  divergence  centered  on  a  modified  cylindrical 
expansion  of  the  pressure.  Essentially,  the  wave  was  allowed  to  expand  at  the  shock  velocity  induced  by 
the  impacting  plate  until  the  shock  wave  at  the  centerline  encountered  the  casewall/explosive  interface. 
A  circle  was  then  fit  to  the  centerline  explosive/cascwall  interface  and  the  radial  (x.y)  position  of  the 
plate.  The  force  was  then  assumed  to  expand  to  the  arc  length  given  in  this  fashion.  This  can  be 
expressed  as: 


Peqv  =  Po*(arclng/di)  (26) 

where  di  is  the  radial  extent  of  the  shock  wave  (equivalent  plate  radius). 

However,  the  interface  pressures  calculated  in  this  fashion  were  high  relative  to  the  hydrocode 
calculations.  A  better  approximation  and  the  one  currently  in  use  is  given  as: 

Peqv  =  Po*(cos  xi/rp)  (27) 

where  rp  is  the  radius  of  the  pressure  pulse.  As  the  impact  becomes  more  one  dimensional,  the  pulse 
radius  becomes  large.  Since  the  plate  diameter  is  limited  to  the  radius  of  the  acceptor  round,  hence  finite, 
the  interface  pressure  approaches  the  impacting  pressure.  The  solution  is  actually  performed  in  two 
iterations.  With  the  pressure  calculated  at  the  interface,  a  new  shuck  velocity  can  be  calculated.  Because 
this  velocity  is  calculated,  a  new  available  time  is  determined.  Figure  19  is  a  plot  of  an  AFX-1 100 
loaded  MK-82  in  a  side-by-side  impact  scenario  at  various  expansions.  The  pressures  plotted  were 
calculated  m  the  explosive  using  both  the  OTl*HULL  code  and  Equation  27. 
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Figure  19.  Plol  of  AFX-1 100  Sidc-By-Sidc  Impact  Scenario  at  Various  Expansions 

Included  in  the  model  is  the  ability  to  insert  a  liner  in  the  cascwall.  Again,  19  shows  the 
results  of  including  the  liner.  It  is  interesting  to  note  that  without  the  liner,  donor  cascwall  expansions 
greater  than  3.5  cm  produced  detonations  in  side-by-side  impact  calculations  involving  AFX-1  lOO 
explosive.  Tests  of  this  same  event  produced  no  detonation  in  side-by-side  configurations.  With  the 
addition  of  the  2.5  mm  thick  liner,  the  model  recreated  the  test  results.  In  addition,  the  pressures  so 
caleulated  were  almost  identical  to  previous  coarsely  zoned  calculations  without  a  liner  that  r.iso 
reproduced  the  test  results.  For  small  expansions,  the  isentjope  pressure  may  be  sii’nincara.  This 
pressure  is  added  such  that  the  total  interface  pressure  is  given  as; 

P  =  Pcqv  +  Pjwl.  (28) 

The  assumption  is  that  the  iscniropc  pressure  pulse  duration  at  the  peak  is  of  the  order  of  the  reaction 
zone  length,  which  is  of  almost  zero  duration  relative  to  the  ease  pulse.  The  pressure  pulse  is  then 
triangular  with  a  duration  at  the  base  given  by  Equation  23. 

The  hydrocodc  results  indicate  that  for  the  pulse  controlled  by  the  plate  thickness,  a  one-dimensional 
region  in  the  explosive  docs  exist  primarily  due  to  the  shock  velocity  difference  between  the  undclonated 
explosive  and  the  metal  casewall.  After  this  region,  the  expanding  wave  equation  by  Green  provides 
good  results. 

With  the  peak  pressure  and  duration  in  the  explosive  determined,  the  likelihood  of  sympathetic 
detonation  can  be  determined  in  many  ways.  The  model  currently  uses  the  critical  diameter  method  to 
evaluate  sympathetic  detonation.  If  the  minimum  go  pressure  recorded  in  the  wedge  test  extends  beyond 
the  critical  diameter  of  the  wedge  test,  a  go  is  assumed. 
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SECTION  III 


DEVELOPMENT  OF  THE  EXPLOSIVELY  DRIVEN  FLYER  PLATE  TEST 

Information  provided  in  Reference  6  was  u:ed  to  devise  a  method  for  accelerating  and  launching  flyer 
plates  into  explosive  targets.  Through  hydrocode  analysis  and  test  constraints,  the  most  feasible 
configuration  consisted  of  a  plate  180  mm  in  diameter  with  a  thickness  up  to  12.7  mm.  In  order  to 
recreate  the  boundary  conditions  observed  in  the  pallet  test,  the  velocity  range  would  be  between  1.4  and 
2.2  km/sec. 

The  basic  requirement  consists  of  using  thick  walled  cylinders  greater  than  76  mm  thick  to  ensure 
shock  wave  rarefactions  from  the  free  surface  do  not  interfere  with  the  plate  acceleration.  In  addition,  the 
pressure  distribution  across  the  plate  face  must  be  as  uniform  as  possible  to  prevent  distortion.  During 
the  detonation  pro;;css  the  center  of  the  explosive  is  the  highest  pressure  region.  Therefore  the  head 
height  of  the  explosive  is  important  in  determining  the  flatness  of  the  plate.  If  the  length  of  the  charge  is 
loo  short,  the  detonation  wave  will  be  very  curved.  If  this  wave  impacts  the  plate,  the  plate  will  be  bent 
accordingly.  Another  technique  that  is  used  to  produce  a  flat  plate  is  the  use  of  an  air  gap  between  the 
explosive  and  tht:  flyer  plate.  If  the  detonation  wave  is  very  curved  when  it  moves  into  the  air  gap,  it 
immediately  slows  down  but  the  edges  of  the  detonation  wave  are  still  traveling  at  the  detonation  velocity 
of  the  explosive.  Therefore  they  catch  up  to  the  shock  front,  and  the  wave  that  impacts  the  plate  is  now 
flat.  Figure  20  is  a  hydrocode  calculation  demonstrating  what  happens  to  a  flyer  plate  when  these 
techniques  are  ivat  applied.  For  this  particular  calculation  the  casewall  of  the  cylinder  is  very  thin,  13 
mm,  and  there  i>  no  air  gap  between  the  plate  and  the  explosive.  Notice  the  curvature  of  the  flyer  plate. 


Figure  20.Hydrocode  Calculation  of  Thin-Wallcd  Cylinder 
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Hydrocodcs  were  used  extensively  in  the  development  of  the  flyer  plate  test.  They  helped  to  identify 
all  aspects  of  the  dctonalion/shock  wave  interaction  and  allowed  a  look  at  different  combinations  of  air 
gap,  head  height,  and  flyer  plate  thickness. 

Figure  21  is  a  hydrocodc  calculation  showing  the  benefits  of  having  a  thick  cylinder  50  mm  wall 
thickness  and  an  air  gap  between  the  explosive  and  the  (Iyer  plate.  Notice  the  flatness  of  the  flyer  plate. 


I'igure  21.  Hydrocode  Calculation  of  Thick -Walled  (Cylinder 
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The  explosive  lhal  was  used  to  launch  the  flyer  plates  was  AFX- 1 100;  it  consists  of  TNT,  wax,  Al 
(66/16/18  percent  by  weight).  It  was  picked  because  it  had  a  low  detonation  pressure  123  kbars,  which 
would  help  prevent  spall  from  occurring  from  the  flyer  plate.  Pentolite  was  used  as  a  booster  to  avoid 
overdriving  the  AFX-1 100  main  charge.  The  detonation  train  consists  of  a  50  x  25  mm  pentolite  booster, 
25  X  25  mm  A-5  pellet,  and  an  RP-83  detonator. 

For  the  test  setup  the  flyer  plate  was  epoxyed  into  a  250  mm  long  by  200  mm  diameter  steel  cylinder 
lhal  had  a  casewall  thickness  of  13  mm.  The  air  gap  device  is  placed  first  behind  the  flyer  plate  and  then 
the  explosive.  With  the  small  can  loaded,  the  next  step  is  to  place  it  into  a  larger  (400  mm  long  by  300 
mm  outside  diameter  with  a  SO  mm  wall  thickness)  steel  cylinder.  By  having  the  smaller  steel  cylinder 
already  loaded,  it  greatly  simplifies  and  prevents  alignment  error  from  occurring  during  the  assembly  al 
the  range.  The  face  of  the  flyer  plate  in  the  small  cylinder  should  be  placed  76  mm  from  the  end  of  the 
bigger  steel  cylinder.  This  distance,76  mm, was  a  comprise;  it  was  a  trade-off  between  having  the  plate  too 
far  back  into  the  barrel,  which  would  cause  damage  to  the  plate  from  the  product  gases, or  totfar  out  with 
no  conoibulion  from  the  product  gases  at  all.  Figure  22  ^ows  the  system  now  being  used.  The  flyer 
plate,  explosive,  air  gap  device  and  steel  cylinders  have  all  been  precision  machined  to  +/-  0.127  mm  for 
a  form  fit. 


^ -  400  mm  - ^ 


Figure  22.  Energetic  Materials  Branch  Flyer  Plate  System 


The  thickness  of  the  flyer  plate  has  been  found  to  be  a  very  critical  element  in  determining  the 
pressure  induced  inside  the  acceptor  charge.  The  proper  thickness  of  the  flyer  plate  was  determined  by 
hydrocode  calculations  and  measurement  of  post-test  bomb  fragments.  Both  the  hydrrKrodcs  and  the 
post-test  data  verified  lhal  the  fragments  had  thinned  down  from  10  to  6  mm.  Therefore  the  flyer  plate 
thickness  that  will  be  used  for  this  investigation  will  be  6  mm.  The  main  method  for  acquiring  plate 
velocity  and  integrity  data  is  flash  radiography.  This  method  is  used  to  look  inside  the  fiicball  through 
the  products  and  x-ray  the  plate  in  flight.  The  plate  is  going  very  fast  (2  km/scc)  ,  so  the  flash  x-ray 
works  well  in  this  lime  frame.  It  verifies  lhal  the  plate  is  flying  flat  with  no  spall.  The  x-ray  system  used 
will  be  the  450  kev  system.  Two  x-ray  heads  will  be  used  so  that  the  velocity  can  be  determined.  The  x- 
ray  film  and  the  heads  are  very  close  to  the  exploding  cylinder,  and  therefore  protection  of  the  heads  and 
a  film  packet  that  can  absorb  air  shcKk  are  important.  For  all  of  the  tests,  the  distance  of  100  mm  from 
the  end  of  the  large  steel  cylinder  was  the  impact  point.  This  was  where  the  plate  no  longer  had  the 
expansion  products  conuibuiing  to  it  and  the  plate  is  flying  flat.  Tlie  first  x-ray  head  is  positioned  at  the 
100  mm  distance.  The  test  setup  is  shown  in  Figures  23  and  24. 
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Figure  24.  Side  View  of  Uic  Flyer  Plate  Test  Setup 

To  reduce  the  overall  plate  velocity,  the  length  of  the  explosive  used  for  launching  was  reduced. 

With  the  reduction  in  length  of  this  charge,  a  large  air  gap  was  used  to  help  reduce  the  effects  from  the 
very  curved  detonation  wave.  However,  as  the  air  gap  was  increased,  the  shock  wave  in  the  steel 
cylinder  was  given  enough  time  to  reach  the  flyer  plate.  The  shock  wave  produced  a  high  pressure  region 
around  the  edges  of '  'le  flyer  plate,  which  caused  it  to  take  the  shape  of  a  saucer.  Figure  25  is  a 
hydrocode  calculation  verifying  this.  Figure  26  is  an  actual  flash  radiography  of  the  plate.  The  solution, 
as  shown  in  Figures  27  and  28,  was  to  isolate  the  plate  from  the  cylinder  wall  using  a  plexiglass  spacer. 
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Figure  25.  Hydrocode  Calculation  Predicting  High  Pressure  Region  on  Edges  of  Plate 


Figure  26.  Flash  X-Ray  of  the  Flyer  Plate  Showing  High  Pressure  on  Edges  of  Plate 
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Figure  27.  Edge  Effects  Elimination  Device 


Figure  28  shows  the  hydrocode  calculation  with  the  steel  flyer  plate  no  longer  in  contact  with  the  steel 
casewall.  Notice  the  flatness  of  the  flyer  plate. 


Figure  28.  Hull  Calculation  With  Flyer  Plate  Not  in  Contact  with  Steel  Cylinder 

With  this  data  a  new  flyer  system  was  developed  as  shown  in  Figure  29.  The  flash  x-ray  data  is 
shown  in  Figure  30.  This  flyer  plate  is  going  2  km/sec. 
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Figure  29.  New  Flyer  Plate  Launch  System  for  Low  Velocity  Plates 


Figure  30.  Flash  X-Ray  of  Flyer  Plate  With  Unifonn  Pressure  Distiibution  Across 
the  Back  Surfa  of  the  Plate 

Having  solved  the  problem  of  reducing  flyer  plate  deformation  at  various  heights,  a  scries  of  tests  was 
conducted  to  relate  explosive  height  to  plate  velocity.  Hydrocodes  previously  conducted  had  shown  the 
region  of  interest  to  be  1-2  km/sec  (Figure  31).  A  calibration  curve  consisting  of  for  shots  is  shown  in 


Figure  31.  AFX  644 1  Calculation  of  the  Detonation  of  the  Donor 
Bomb  With  a  Flyer  Plate  Foiming  at  75  ps 
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Flyer  Velocity  -vs-  E)?)losive  1hid<ness 
28  Eqn  196  y=((a+cxMl+bx))  r2=0.995263479 
0=0.0086121857  b=0010613488 
c=00351 27037 


Figure  32.  Plot  of  Flyer  VelocUy  versus  Explosive  Thickness 
The  equation  that  was  derived  from  this  plot  is  shown  below. 

y  =  ((a  +  cx)/(l  +  bx)) 

where  0.995,  a  =  0.0086,  b  =  0.0106,  c  =  0.03512. 


(29) 
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SECTION  IV 


SHOCK  TO  DETONATION  TRANSITION  IN  CYLINDRICAL  CHARGES 

Wi(h  the  flyer  plate  velocity  curve  established,  live  acceptors  could  now  be  tested  using  the  flyer  test. 
The  live  acceptors  were  explosively  filled  right  circular  cylinders  with  piezoelectric  pins  embedded  into 
the  explosive.  The  cylinders  are  1018  cold-roUed  steel  250  mm  long  by  200  mm  outside  diameter  and 
180  mm  inside  diameter.  The  acceptor  cans  have  17  pin  holes  drilled  along  the  long  axis  and  centered  on 
the  can.  Piezoelectric  pins  are  placed  in  these  holes  in  a  stair-stepped  fashion  away  from  the  inside 
surface  of  the  acceptor  casewall.  A  gauge  is  used  to  obtain  proper  spacing  of  the  pins  from  the  casewall. 
The  17  pins  make  up  two  channels  of  data  with  the  center  pin  touching  the  inside  surface  of  the  acceptor 
casewall.  A  drawing  of  the  setup  is  shown  in  Figure  33.  Figure  34  is  a  side  view  of  the  test  setup;  notice 
the  piezoelectric  pins  embedded  in  the  acceptor  cylinder. 


Acceptor  Cylinder 


Flyer  Plate 


Figure  3-3.  Flyer  /Acceptor  Setup  With  Pins  in  Place 


Piexoeklric  Piai 


Explosively  Filled  Cylinder 

i 
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Figure  34.  Side  View  of  the  Flyer/Acceplor  Setup 
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A  second  method  of  obtaining  data  was  the  use  of  a  250  x  250  x  25  mm  thick  witness  plate  under  the 
acceptor  charge.  The  armor  shows  if  the  acceptor  cylinder  detonated.  If  the  acceptor  does  not  detonate, 
the  armor  will  be  unmarked.  However,  this  does  not  mean  the  acceptor  cylinder  did  not  undergo  some 
type  of  reaction.  The  other  instrumentation  gives  this  kind  of  detaU.  Figure  35  is  a  photograph  of  the 
overall  test  setup. 


Figure  35.  Overall  Test  Setup 
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SECTION  V 


EXPERIMENTAL  RESULTS 

Figure  36  dcpicis  the  comparison  of  the  simplincd  model  data  to  the  experimental  method.  The  test 
series  utilized  an  insensitive  high  explosive  (IHE)  AFX -6441  in  the  acceptor  charges.  The  experimental 
data  was  gathered  using  piczociccuic  pins.  All  of  the  data  was  ruicd.  The  IHE  used  in  the  acceptor 
exhibited  a  go-no-go  point  ui  a  pressure  of  approximately  70  kbars^which  corresponds  to  a  shock  velocity 
of  4.04  mm/iiS| which  results  in  a  detonation.  The  run-up  equation  derived  from  ELSGT  tests  is  log  X*  = 
•  3.24  log  P*  -f  S.2I .  A  shuck  velocity  of  4.2  mm/fis  corresponds  to  a  pressure  of  81  kbars,which  results 
in  a  run  of  185  mm.  The  detonation  velocity  of  this  explosive  is  7.1  mm/ps.  From  Figure  36, by  120  mm 
the  transition  to  detonation  is  almost  complete.  The  implication  is  that  the  induced  pressure  is  actually 
higher  or  the  pop  plot  is  in  error.  Tlic  model  predicts  a  detonation  for  the  1 .9  mm/ps  flyer  against  this 
explosive  using  a  12.7  mm  critical  diameter  and  no  detonation  at  1.48  mm/ps.  Interface  shock  velocities 
predicted  by  the  model  arc  upproximaicly  4  percent  higher  in  all  cases. 


Figure  36.  Comparison  of  the  Simplified  Model  With  the  Experimental  Data 
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Figure  37  is  ihe  manganin  stress  dnia  for  the  1.9  flyer  velocity.  The  prediclcd  pressure  at  this 

point  (13  mm  into  the  explosive)  was  55  kbars.  The  constantan  strain  gauge  revealed  that  the  manganin 
gauge  started  receiving  a  non-planar  stress  wave  approximately  l(X)  nanoseconds  into  the  event.  It  is 
believed  that  the  true  peak  pressure  signal  was  lost  due  to  this  problem.  The  technique  Is  still  being 
perfected,  and  more  tests  will  be  performed  to  obtain  the  actual  pressure. 

The  pulse  duration  for  this  test  as  shown  in  Figure  37  is  given  as  0,986  ps.  This  compares  TOA 
predicted  pulse  duration  of  1.07  ps.  It  is  not  known  the  exact  effect  the  non-planer  wave  has  on  pulse 
duration. 

MANGANIN  STRESS  GUAGE 


Figure  37.  Manganin  Gauge  Data 
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SECTION  VI 


DISCUSSION  AND  CONCLUSIONS 

The  purpose  of  ihis  research  effort  was  to  develop  a  means  of  predicting  the  likelihood  of  sympathetic 
detonation  in  the  storage  configuration  for  explosive  fills  under  development.  Since  it  is  not  clear  which 
mechanism  or  combination  of  mechanisms  is  responsible  for  sympathetic  detonation  (these  may  vary 
from  explosive  to  explosive),  the  goal  was  to  develop  a  test  that  accurately  reproduces  the  conditions 
observed  in  the  storage  configuration.  If  the  boundary  conditions  can  be  reproduced,  then  those 
mechanisms  should  exhibit  themselves.  While  it  is  possible  to  instrument  rounds  in  the  storage 
configuration,  this  experiment  allows  much  greater  control  of  the  acceptor  environment. 

Results  of  the  model  indicate  interface  pressures  approximately  12  percent  above  those  recorded  in 
the  experiment.  However,  at  12.7  mm,  the  calculated  pressure  for  the  1.9  mm/ps  impact  test  is  83  kbars 
or  2  percent  above  the  experimental  result.  Assuming  the  nmup  was  induced  by  the  calculated  interface 
pressure  of  93  kbar,  the  predicted  runup  is  1 18  mm,  significantly  closer  to  that  observed  in  the 
experiment  than  would  be  predicted  by  an  81  kbar  initiation.  The  model,  as  mentioned  previously, 
predicts  sufficient  stimulus  to  achieve  detonation.  Only  the  physical  size  of  the  acceptor  charge  prevents 
this.  This  suggests  the  interface  values  predicted  by  the  model  are  reasonably  accurate  but  the  rate  of 
shock  wave  decay  may  be  too  low.  The  calculated  pressures  at  12.7  mm  into  the  explosive  for  the  1.48 
mm/|is  impact  point  is  S 1  kbar.  The  one-dimensional  region  in  the  explosive  varies  but  was  generally 
controlled  by  lateral  relief. 

Differences  between  the  model  and  the  experiment  may  be  due  to  curvature  in  the  impacting  plate, 
small  errors  iti  impact  velocity,  or  differences  in  the  Hugoniot  modeled  and  that  tested.  Regardless,  the 
model  predicts  with  a  high  degree  of  accuracy  the  interface  condition  in  the  round-to-round  acceptor 
experiment,  giving  the  explosives  formulator  a  quick  and  accessible  tool  to  determine  the  likelihood  of 
sympathetic  detonation.  Similarly,  the  flyer  plate  test  has  proven  u>  he  a  versatile  and  valuable  tool  in 
assessing  safety  margins  or  lack  thereof  of  rounds  in  the  storage  configuration. 
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To  whom  it  may  concern, 

Greetings,  my  name  is  J.  Gregory  Glenn,  you  have  recently  received  a  report  titled  Simulating 
Sympathetic  Detonation  Effects  (Report  #  WL-TR -93-7030).  In  my  hurry  to  finish  the  report  I  left  out  a 
very  important  part.  The  part  that  I  am  referring  to  is  the  acknowledgments  located  on  page  iii.  1  am 
sending  the  proper  acknowledgments  and  am  asking  you  to  please  replace  this  with  the  one  in  the  report. 
I  would  fill  remiss  if  I  did  not  acknowledge  the  help  that  I  received  from  these  people.  Thank  you  for 
your  help  in  this  matter. 
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